Dynamics of kicked particles in a double-barrier structure 
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We study the classical and quantum dynamics of periodically kicked particles placed initially 
within an open, double barrier structure. This system does not obey the Kolmogorov-Arnold-Moser 
(KAM) theorem and displays chaotic dynamics. The phase space features induced by non-KAM 
nature of the system leads to dynamical features such as the non-equilibrium steady state, classically 
induced saturation of energy growth and momentum filtering. We also comment on the experimental 
feasibility of this system as well as its relevance in the context of current interest in classically induced 
localisation and chaotic ratchets. 
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I. INTRODUCTION 

Periodically kicked rotor is a popular model that has 
served as a paradigm to understand Hamiltonian chaos 
both in the classical and quantum regime [l]. This 
was originally introduced as a simple model for dynam- 
ical chaos but was sufficiently general enough to cover 
many physical situations. For instance, problems like 
the Hydrogen atom in microwave fields and motion of a 
comet around the sun driven by a suitable planet can 
be reduced to that of kicked rotor Q. This system is 
also paradigmatic for another important reason; it obeys 
KAM (Kolmogorov- Arnold-Moser) theorem j3[. This 
implies that, as a control parameter is varied, the tran- 
sition from regularity to chaos occurs progressively by 
breaking of invariant curves in phase space. Once all 
the invariant curves are broken down, diffusive global 
transport of particles in phase space becomes possible. 
In the corresponding quantum regime, this classical dif- 
fusive transport is inhibited by the onset of dynamical 
localisation !4J. This was experimentally realised in the 
laboratory with cold atoms in optical lattices [5| and is 
the basis for theoretical and experimental realization of 
chaotic ratchets in recent times [6|. 

On the other hand, there are other physical systems 
that exhibit classical chaos but violate the KAM theo- 
rem, the so-called non-KAM systems. This class includes 
the kicked harmonic oscillator (KHO) Q and the kicked 
particle in infinite square well potential [1, ■ In both 
these cases, when a parameter is varied, the invariant 
curves are replaced by stochastic webs [7|, an intricate 
chain of islands and globally connected channels, through 
which particle transport becomes possible. Non-KAM 
type is also relevant for an important class of physical 
systems, namely the dynamics of particles in quantum 
wells and barrier structures. Till date, non-KAM sys- 
tems have been experimentally realised in semiconduc- 
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tor superlattices in tilted magnetic fields in which the 
enhanced conductivity could be attributed to non-KAM 
chaos [111 ] . Further, measurement of Lochschmidt echo 
using a non-KAM system, namely, the ion trap with har- 
monic potential in the presence of a kicking field [HJ has 
also been proposed. Inspite of this, very few non-KAM 
systems have been investigated and they have not been 
explored in sufficient details. 

Another motivation for this work stems from the con- 
siderable interest in recent times in the dynamics of con- 
densates placed in finite box-type potentials acted upon 
by a periodically kicking field. In a recent experiment, 
Henderson et. al. [H[ have constructed a quasi- ID finite 
box, using a combination of optical and magnetic trap, 
with the Bose-Einstein condensates in the box receiving 
periodic kicks. This set-up was used to study the effect of 
atomic interactions on the transport of BECs. In place 
of the dynamical localisation they observed a classical 
saturation in the energy of BECs due to a balance be- 
tween the energy gained from kicks and the energy lost 
by leakage of BECs over the finite barrier Apart 
from this, a series of experiments [l4| that studied the 
transport of BECs in the presence of disordered potential 
have reported such classically induced energy saturation 
effects. Then the question is if it is possible to observe 
such classically induced energy saturation in chaotic sys- 
tems without inter-particle interactions and what would 
its mechanism be ? We show that kicked particle in finite 
well type potential that we study in this paper shows this 
feature and we discuss its mechanism. It is also relevant 
to point out that following the achievement of BECs in 
the optical box trap , theoretical investigations of res- 
onance and anti-resonance behaviour and its relation to 
the KAM and non-KAM type dynamics for BECs in ID 
infinite well have also been performed [l6| . Further, ex- 
periments exploring the interface of nonlinear dynamics 
of electrons in ID quantum well irradiated at terahertz 
frequencies have already been reported [l7j . 

Although non-KAM type dynamics is a generic fea- 
ture in physical systems such as the potential wells not 
much work has been done on this class of problems. How- 
ever, on the theoretical front, infinite square well po- 
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FIG. 1: Schematic of the stationary part of the potential. The 
width of each barrier is b. The well region has width 2a. 

tential confining a delta-kicked particle has been studied 
d, [l(| ■ In the light of recent attempts to study chaotic 
ratchets @, in both the classical and quantum sense, it 
would be useful to open the potential to allow for particle 
transport. This could lead to chaotic ratchets that can 
utilise non-KAM features of its classical system for di- 
rected transport. In this paper, we study a periodically 
kicked particle initially held in between a finite double 
barrier structure. Double barrier heterostructures play 
an important role in electronic devices that use resonant 
tunnelling diodes [lH though without the kicking poten- 
tial. Primarily we present numerical explorations of this 
problem to study its rich dynamical features. In the next 
section, we introduce our model and in subsequent sec- 
tions we discuss the classical and quantum dynamics of 
this system. 



with E c = m\ 2 /2T 2 leads to a new dimensionless Hamil- 
tonian 

2 00 

H=y + V sq + ecos(x) *(*-»)• (3) 

n— — oo 

In this, V sq = Vo [Q(x -<j) + b + B,n)-<d(x-(j) + Rtt) 
+0{x - 4> - Rtt) - Q(x - <f> - Rtt - b)} with R = 2a/ X 
being the ratio of the distance between the barriers to 
the wave length of the kicking field. The classical dynam- 
ics of the system depends upon five parameters, namely, 
e, R, fo, Vq and <j). Of these, R, b and <p determine the po- 
sitions of discontinuities in the potential (position of the 
wall boundaries) at B = {—xi — b, —xi, x r , x r + b} where 
xi = —Rtt + (j> and x r = Rtt + cf>. Note that if <f) = 0, 
then X\ = —x w and x r = x w with x w — Rtt. Thus, 
the qualitative nature of the classical dynamics depends 
on the positions of wall boundaries collectively denoted 
by B, the kick strength e and potential height Vq. In 
this paper (except in section 111(C)), we set <\> = which 
makes the potential symmetric about x = 0. It is useful 
to write Eq. [3] as 

H = H + V sq {x), (4) 

where H = ^- + e cos (x) J2^L-od ~ n ) leads to stan- 
dard map defined on the infinite plane. Note that if 
V sq (x) = Vq, a constant, then the Hamilton's equations 
will not have the potential term and the dynamics would 
be completely governed by Hq. 

III. CLASSICAL DYNAMICS 



II. THE KICKED PARTICLE IN DOUBLE 
BARRIER 

We consider the dynamics of a non-interacting par- 
ticle initially located in between two potential barriers 
each of height Vo and width b distant 2a apart (see Fig. 
[1]). The particle is further subjected to flashing £-kicks 
of period T generated by a spatially periodic potential 
field of wavelength A. Amplitude e of the kicking field 
is generally referred to as kick strength. The classical 
Hamiltonian of the system is 



^+0) J2 s(t-nn (1) 



~ p 2 ~ _ 

H = 1- V sq + e cos 

2m V A 



where V sq = Vq 



-Q(x- 



Q(x + b + a)- Q(x + a) + &(x - a) 
a — b) , 0(.) is the unit step function and <p is 



the phase of the kicking field. The set of canonical trans- 
formations given by, 
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A. Classical map 

The Hamiltonian in Eq. Q is classically integrable 
for e = 0. This corresponds to free motion in the pres- 
ence of two potential barriers and hence it is possible to 
obtain a transformation to action-angle variables. For 
e > 0, the system is non-integrable and can even display 
abrupt transition to chaotic dynamics with mixed phase 
space depending on the values of R, b and <j). It is conve- 
nient to think of the system as being entirely governed by 
Hq and then incorporate effect of discontinuities in V sq 
through appropriate boundary conditions. This leads to 
the following map, 



Pn = Pn-i + esin(a;„_i), 



n 



(5a) 



(5b) 



Equation [Sa] represents the effect of Ho and is iden- 
tical to the standard map. In Eq. I5bl the operator 
7Z = TZk ■ ■ ■ 72-2 7£i represents the effect due to k encoun- 
ters of the particle, in between two kicks, with the discon- 
tinuities of V sq at positions represented by B\ , Bi , . . . 
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respectively. Depending on the energy, each of these 
k encounters could either be a reflection (sign of mo- 
mentum changes) or refraction (magnitude of momentum 
changes) at £>,; € B, i = 1, 2, ....k. 

The map in Eq. [5] would be complete if the oper- 
ator IZi, that incorporates effect of i th discontinuity en- 
countered, is explicitly written down. Between successive 
kicks applied at times n and n + 1, we denote the state 
of the particle after incorporating effect of ith encounter 

with a boundary Bi by ( ^" ^ . We define ] x\ , a4 [ with 

i = 0,1... k as the path, starting from x % s , a particle 
would traverse between the two kicks after encountering 
i th discontinuity if there were no discontinuities to be 
faced till the next kick. For i = 0, x l s would simply be 
x n -i and would be equal to Bi for i > 0. x° and p° 
to be used in boundary conditions would simply be x n 
and p n obtained directly from Eq. I5al Boundary condi- 
tions defined by Eq. [B] below are applied k times until 
]xg,x£[ n B = 0. If E n denotes the energy of the sys- 
tem at nth kick, then for E n < Vq (reflective boundary 
condition), we obtain 

^ / a:*" 1 \ (2B i -xt l \ 

Kil = . (6a) 

{Pn 1 ) V -Pn' J 

For E n > Vq (refractive boundary condition), we get 
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FIG. 2: (Color online) Stroboscopic Poincare section (black) 
for R = 0.95, e = 0.15, Vo = 0.5, cj> = and b = 0.5. All the 
continuous curves (in color) marked C\ to C@ are for the cor- 
responding standard map with kick strength 0.15. The black 
box at position x = ±x w indicates the width b of the barrier. 
The solid circles (in red) show a trajectory starting from A\ 
until it exits the potential well at Ag. The time ordered se- 
quence of the trajectory is A\ to A2, reflection at —x w , A3 to 
A4,, reflection at x w , A5 to ^6, cross the boundary at x w , A7 
to As, cross the boundary at x w + b, exit the potential at Ag. 
See text for details. 
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ZVoVdiff 
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(6b) 



In this, we have used Vdiff — V{x l ~ 1 ) — V{x l j T 1 ). Thus, 
the dynamics of system in Eq. [3] can be described by the 
standard map with — 00 < x n ,p n < 00 (Eq. [5J subjected 
to potential barriers (Eq. I5b[) . Notice that by putting 
Vo = in Eq. I6b[ we obtain IZi — I for all i, where I is 
the identity matrix of order 2. Then 1Z = I and, as ex- 
pected, Eq. [5]reduces to standard map for Vo = 0. Thus, 
the transformation ([6]) can be viewed as deviation from 
standard map dynamics induced after each encounter of 
the particle with the a discontinuity of potential V sq . 



B. Phase space features 

Figure [2] shows a stroboscopic section obtained by 
evolving the map in Eq. [5] for uniformly distributed ini- 
tial conditions in x G ( —PcPc), where 
p c = \/2mVo is the minimum momentum required for 
barrier crossing. In this paper, we have chosen kick 
strength e << 1 such that the corresponding standard 
map displays only KAM curves. Firstly, a striking fea- 
ture is the absence of invariant curves and the appearance 



of a mixed phase space. This is in stark contrast with 
the standard map which displays mostly quasi-periodic 
orbits for kick strengths of this order. This figure also 
shows snap shots (solid circles in red) of trajectory in- 
between successive encounters with the discontinuities at 
B. Clearly, the evolution between two successive encoun- 
ters with the boundaries is confined to a trajectory that 
is identical with one of the quasiperiodic orbits of the 
corresponding standard map (obtained from Eq. [S]with 
Vo = 0) shown as continuous lines in the figure. Due 
to V sq , particle breaks away from one quasiperiodic orbit 
and joins another at each encounter with the boundaries. 
This leads to absence of quasiperiodic orbits and develop- 
ment of mixed phase space comprising intricate chains of 
islands embedded in chaotic sea. We illustrate the effects 
of discontinuities in Fig. [5] by following a typical initial 
condition marked A\ in the chaotic layer. This evolves to 
A2 on the invariant curve C2 of the corresponding stan- 
dard map. After a long time, this point appears on the 
curve C\ and goes from A3 to A4. After a reflection at 
—x w , it goes from A5 to Aq on C3. Then it shifts to 
the barrier region (x w ,x w + b) and moves on C7 from 
A7 to Ag,. Depending on the winding number of the or- 
bit in (x w , x w + b), the particle could have gone back in 
to region between the barriers or escape from the finite 
well. In the present example, it makes its escape out of 
two barrier structures and its state meets the curve C5 
at Ag. Once the particle has escaped, its state evolves 
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on same curve as n -> oo. Thus, system displays KAM 
behaviour for |x| > x w + b. 

The absence of quasiperiodic orbits can be attributed 
to the non-analyticity of V sq which violates the assump- 
tions of KAM theorem. Thus, the non-KAM nature of 
system leads to onset of chaos even for e < 1. The 
initial conditions starting from chaotic layer will diffuse 
in momentum space. Some of these initial conditions 
which reach the set of quasiperiodic orbits C (//) (// being 
the winding number) of the corresponding standard map 
which overlaps the region \p\ > p c can escape from the 
finite well. As \x increases, this overlap also increases and 
hence the escape probability is larger. This implies that 
there must exist /z c such that the states on any C(p), 
with /j, > fa, will definitely cross the barrier and escape 
from the well. These orbits do not encounter the discon- 
tinuities in the potential multiple times and hence the 
energy of the particles evolving on such quasi-periodic 
orbits will not diffuse. Figure [2] also shows the trajec- 
tory of a particle (open circles in red on the curves C5 
and Ce) in such non-diffusive region. As seen in Fig. [21 
the discontinuities at x w and x w + b relocate the incom- 
ing particle from C5 (/ig) to another orbit Cq (hs), where 
/U5 and fie are their winding numbers, respectively. As 
shown in Appendix B, when b — > 0, the effect of these 
discontinuities decreases and deviation between between 
two orbits measured as {fi& — /U5) — > 0. This results in 
the appearance of regular orbits (see Fig. [3]) identical to 
those of the standard map except that the former have 
imperceptible discontinuities wherever there is a discon- 
tinuity in potential. In other words, refraction becomes 
identity operation as b — > 0. Thus, the system shows 
regular dynamics outside region enclosed between curves 
C± (a*c) (see Fig. [3]) as b — > 0. Note that the limits on 
chaotic phase space in terms of fi on positive and nega- 
tive sides of momentum are identical due to assumption 
that = 0. Limits on the chaotic phase space would ex- 
ist even otherwise, though these would not be identical 
on both sides of p = 0. The discussions in this sub- 
section can be summarized as follows ; we can define a 
phase space region Ai (\x\ < x w + b; \p (x) \ < p (x; He)), 
such that system has mixed phase space inside Ai in 
general and regular dynamics outside it. Here, p (x; /i c ) 
is momentum of any state on the curve C+ (/Lt c ) at posi- 
tion x. In Fig. [3l a close numerical approximation of the 
region Ai is highlighted by the red dashed line. 

We remark that for b — > 0, the phase space structures 
inside At are identical to those of well map that describes 
the dynamics of 5- kicked particle in an infinite well Q. 
This is to be expected since the well map has only re- 
flective boundaries for \p\ < 00. Further, the well map is 
hyperbolic for R < 0.5 for any e > 0. The Hamiltonian in 
Eq. [I] also displays complete chaos for R < 0.5 inside Ai. 
This is seen in Fig. [3] as no regular structures are visible 
in this region to the accuracy of our calculations. The 
region defined by Ai is determined by the positions of 
potential discontinuities B and C±(/i c ). It can be shown 
that C±(fj-c) will remain close to ±p c (— ±\/2mVo) when 
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FIG. 3: Stroboscopic plot (excluding C±(/X{,)) for b — 10~ 3 
and R — 0.5. All the other parameters are same as in Figure 
[2] Dashed line (in red) represents the boundary of region 
Ai. The mild scatter of points just below C+(fj,c) and just 
above C_(/x c ) represent the particles escaping out of the well 
(whose initial states were in Ai). The curves C+(nb) and 
C-(nb) shown here are used in section V(B). 

b — > for any e for which standard map has mostly reg- 
ular phase space. Thus, the extent of chaotic region will 
depend grossly on the positions B and height Vq of the 
barriers only. This implies that it is possible to engineer 
chaos in a desired region by varying these parameters. 

C. KAM-like behavior: Role of symmetries 

In this section, we explore the conditions under which 
KAM or non-KAM type of dynamics can be realized in 
the system. In Eq. [T] the non-analyticity of V sq violates 
the assumptions of the KAM theorem. Hence, generically 
we expect this system to display the signatures of non- 
KAM system such as the stochastic webs instead of quasi- 
periodic orbits and an abrupt transition to chaos. These 
features are shown in Fig. 2]Ja,c,d,f). However, we show 
that even in the presence of non-analyticity in V sq , quasi- 
periodic orbits similar to that in a KAM system can be 
realised if certain symmetry conditions arc satisfied. 

As argued before, until interrupted by the barriers, the 
dynamics is confined to a particular invariant curve of the 
corresponding standard map. We recall that correspond- 
ing to every trajectory C+ of standard map with p n > 0, 
there exists one and only one trajectory C_ with p n < 0, 
such that a particle will evolve on these trajectories in ex- 
actly the same way but in opposite direction. As shown 
in appendix-A, consider the (R, 4>) pairs for which the 
condition 

± Rtt + <j) = I 2tt, I eZ (7) 
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The scaled Planck's constant is K s = =§-?p. This being 
a kicked system, we can obtain the one-period Floquet 
operator, 



U = cxp I — cos x ) exp ( — — 

' n. J V h x 



p T + v„ 



0) 



such that ifj(x,n) = U n t/j(x,0). The classical limit will 
correspond to taking h s — ¥ keeping e = eh s /h constant. 
We calculate the Husimi distribution Q(xq,pq, n) defined 
by 



Q(x ,p Ql n) = \(ip(x, n)\x ,p )\ 2 



(10) 



FIG. 4: Stroboscopic Poincare section for the Hamiltonian 
in Eq. [1] showing the region x £ (— xi,x r ),p £ (— p c ,p c ) for 
b = 0, e = 0.15, Vb = 0.5. The other parameters are (a) 
R = 0.95,0 = (b) R = 1.0,0 = 0, (c) R = 1.05,0 = 0, 
(d) R = 0.45,0 = tt/2, (e) R = 0.5,0 = tt/2 and (f) R = 
0.55,0 = tt/2. 



is satisfied. When Eq. [7J is satisfied, application of TZi 
takes a particle from C+ to C_ and application of 7£j+i 
brings it back to C + . This leads to quasiperiodic be- 
havior in which the particle is confined to a pair of tori. 
This quasiperiodic orbit undergoes smooth deformation, 
just like in a KAM system, until it breaks for large kick 
strengths. Hence we call this KAM-like behaviour for 
its striking resemblance to the qualitative behaviour of a 
KAM system. In general, there exist infinite (i?,0) pairs 
for which KAM-like dynamical behaviour can be recov- 
ered in this system. In Fig@|b,d), we show the sections 
for R = 1, <f> = and R = 0.5, = tt/2 for which KAM- 
like behaviour is obtained. In Fig[4|a,c,d,f), we also show 
cases where Eq. [7] is not satisfied and hence for |p| < p c 
stochastic webs and chaotic regions are seen. 

Symmetry related invariant curves like C+ and C_ are 
due to the symmetry of the kicking field about any x = 
to7t + where m is an integer. It turns out that when 
Eq. [TJis satisfied, kicking field is symmetric about x w and 
X- w . The existence of KAM-like behaviour in presence 
of non-analytic potential can be attributed to existence 
of centres of symmetry of kicking field at — x w and x w . 



IV. QUANTUM DYNAMICS 

In this section, we discuss the quantum simulations 
of the system. We start by writing down the time- 
dependent Schroedinger equation corresponding to the 
scaled Hamiltonian in Eq. ([3]), 



ifr, 



dip 
~dt 



2 8x 2 



sq 



ecosx^^ 8{t 



for a wavepacket at time n. In this we take (x\xq,po) 
as the minimum uncertainty wavepacket. In the semi- 
classical regime, the dynamics in the Husimi representa- 
tion mimics the classical dynamics of the system in phase 
space @. In Fig. El we show the Husimi distribution at 
n = 250 from which one can clearly see that the density 
of Husimi distribution shows pattern similar to classical 
structures shown in Fig. [2j 

Since p and V sq in Eq. [S] do not commute, we first 
divide the duration between successive kicks into N& t 
small time steps and the second term of Eq. (J9j) be- 

n ex p 



comes 



Then, we ap- 
to evolve the system. 



V>. (8) 



h s N At 

ply the split-operator method [it 
We use Fast Fourier transform [2(| to obtain ip (p) from 
if) (x) and vice-verse. In our calculations, we have taken 
N/st ~ 2500, the typical temporal step size is O(10~ 3 ) 
and spatial step size is O(10 -4 ) to ensure that the evolved 
wavepackets converged to at least 8 decimal places. 

The initial wavepacket at n — is located in between 
the two barriers. We choose parameters b and H s for 
which the Husimi distribution (shown in Fig. [5]) closely 
resembles the classical phase space and shows that the 
probability density associated with the initial wavepacket 
will ultimately leave the barrier region by predominantly 
following the classical path rather than by tunnelling. 
Thus, the system stays in the semiclassical regime and 
tunnelling is suppressed. Quite clearly, for such a choice 
of parameters in the semiclassical regime, the classical 
dynamical features would be reflected in the quantum 
dynamics as well. 

In the next two sections, we discuss some interest- 
ing dynamical features, namely (i) the non-equilibrium 
steady state (ii) classically induced suppression of diffu- 
sion and (iii) momentum filtering which primarily arise 
due to co-existence of diffusive (chaotic region (/i < /i c )) 
region and non-diffusive region (regular region (// > /i c )) 
in same non-KAM system. 
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FIG. 5: (Color online) (Top) Husimi distribution for evolved 
wave packet. Initial wave function corresponds to Q(xo,po, n) 
sharply localised inside the chaotic region around (0,0). In 
the grey scale version, the regions with larger values of Husimi 
distribution function are grossly represented by the darker 
areas. It shows that the function decays very steeply outside 
[x—u,,x w ] and acquires negligible values compared to those for 
region inside [x- w , x w \. We have taken h s = 0.0025, R = 0.85, 
b = 0.2, e = 0.15, Vo = 0.5, = 0. (Bottom) Enlarged and 
better resolved view of inset from figure on the top shows path 
followed by probability density outside the barrier region. 



V. DYNAMICAL FEATURES 

A. Non-equilibrium steady state 

In this section, we show that the system in Eq. [1] 
can support non-equilibrium steady state (NESS) for in- 
termediate time scales. We start with initial conditions 
uniformly distributed on a thin rectangular band around 
p = stretched across the well region in between the 
potential barriers. As the kicking field begins to impart 
energy to the system, the particles which absorb suffi- 
cient energy escape from the well. At any time n, the 
mean energy (E)i n of the particles lying inside the well 

is ^^f^i where (.) represents average at time n over the 

classical states (evolved from initial states over n kicking 
cycles) for which —xi < x < x T . In the corresponding 
quantum regime, we have, 



\E)i n — 



ip*(x,n) —ip(x,n) dx (11) 



The effect of the operator p 2 on ip(x, n) can be calcu- 
lated using fast fourier transform and is equal to inverse 
fourier transform of p 2 ip(p,n). Figure |5] shows that ini- 
tially (E}i n increases and after a time scale t r , (E)i n 



saturates to a constant. During this time scale, the be- 
haviour is similar to the classical diffusive regime of the 
standard map. 

The existence of steady state can be understood as 
follows. For the parameters used in Figure [6] the phase 
space in region M. is fully chaotic. As kicks begin to act, 
any localized classical distribution po(x,p) is quickly dis- 
persed throughout this region. The total energy E n of 
the particles in the well region increases. Simultaneously, 
the particles with |p| > p c leave the finite well leading to 
loss of energy. Soon the loss process becomes significant 
and at every kick cycle the energy lost (due to barrier 
crossings) is more than the energy gained from the kick- 
ing potential. Thus, E n begins to decrease. However, 
after the time scale t r , the net energy change and the 
number of particles vary in such a manner as to maintain 
the mean energy (£%„ a constant (apart from fluctua- 
tions). This arises because the normalised momentum 
distribution remains nearly invariant with time as shown 
in Fig. [7] The chaotic mixing inside the well ensures 
that, despite the loss of energetic particles, momentum 
distribution remain invariant. Thus, chaos between the 
barriers is essential to support the NESS. One of the fac- 
tors that determine t r is the rate at which any initial 
distribution of states diffuses in the chaotic region and 
steady state distribution shown in Fig. [7] is achieved. 
This rate increases with e in general. For the present 
case with complete chaos, one expects this rate to be 
proportional to 1/e 2 , just like in the diffusive regime of 
standard map and hence one expects t r oc 1/e 2 . Numer- 
ical results shown in Fig. [6] show a good agreement with 
this gross estimate for t r . 

This steady state holds good until nearly all the parti- 
cles have escaped out and only a fraction q < < 1 remains 
in the well. Based on rate of diffusion in chaotic region, 
we can estimate the time at which this happens to be 
t s cx l/q 2 e 2 . Since q « 1, we get t s >> 1. In the 
semiclassical regime, this mechanism carries over to the 
quantum dynamics as well. Notice that t s is larger than 
other relevant time scales, i.e, t s >> t r > T. Further, t s 
is typically about few hundreds of kick cycles and hence 
we expect this to be experimentally accessible time scale 
as well. On a much longer time scale as t — > oo, all the 
energetic particles escape and the steady state decays 
out. 

Indeed, a similar non-equilibrium steady state has 
been experimentally observed with periodically kicked 
Bose-Einstein condensate in a finite box for strong kick 
strengths ||. These steady states have a classical expla- 
nation. Typically, the standard kicked rotor exhibits en- 
ergy saturation and steady state, for large kick strengths, 
in the quantum regime due to destructive quantum inter- 
ferences. We emphasise that the energy saturation, 
in our model as well as for the BEC in finite box [5j, is 
induced by the classical effects and leaves a trail in the 
semiclassical regime. 

Fig [6] shows that the quantum mean energy (E) in fol- 
lows the classical curve quite closely. These results cor- 
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FIG. 6: Non-equilibrium steady state in the system in Hamil- 
tonian Eq. [TJ The mean energy for the particles held in 
between the double barrier structure. The solid lines are the 
classical results and the symbols correspond to quantum re- 
sults. The other parameters are R = 0.5, b = 0.2, cj> — 0, Vq = 
0.5 and for quantum simulations h s = 0.0025. The solid sym- 
bol (triangle up) marks the time scale t r at which the system 
relaxes to the steady state. 




P 



FIG. 7: (Color Online) Classical steady state momentum dis- 
tribution for e = 0.25 at (a) n = 100 (solid) and (b) n = 200 
(dashed). As seen from Fig. [6jb), the steady state is reached 
at n ~ 30. The other parameters are R = 0.5, Vo = 0.5, b = 
0.2,0 = 0. 



respond to h s = 0.0025 and reflect the behaviour in the 
semiclassical regime. Larger values of e correspond to 
moving away from semiclassical regime towards purely 
quantum regime. Thus, we should expect quantum aver- 
ages to deviate from classical averages in a pronounced 
manner. This is borne out by the numerical results in 
Fig |Bfa,b,c). There is current interest in quantum non- 
equilibrium steady states about which not much has been 
explored until now [2l|. For e >> 1.0, the quasiperiodic 
orbits of the standard map are sufficiently destroyed to 



allow global transport in phase space. Then, particles 
do not have to rely on discontinuities in V sq to diffuse in 
phase space. This leads to unlimited energy absorption 
by the particles between the barriers and NESS is not 
supported. Then, the system essentially works like the 
kicked rotor in the strongly chaotic regime. 



B. Energy saturation and steady state 

As pointed out earlier, the region covered by the quasi- 
periodic orbits C(/i) with \i > fi c is non-diffusive. Now, 
consider curves C± (fib) such that maximum value of \p\ 
for both these curves is equal to p c as shown in Fig. 
[5] Then, all the curves C(fi) with ji < ji^ will have 
p € [— PaPc]- So any state evolving on one of them will 
always get reflected at the barriers. Thus, to escape from 
the finite well, every phase space point in the chaotic re- 
gion must first reach any C(/z) with ji^ < /j, < fi c . As 
time n — ¥ oo, all the particles would have escaped from 
the well and get locked on to one of the invariant curves 
C(/Lt) of the corresponding standard map. Thus, the mo- 
menta of escaping particles settle to a stationary distri- 
bution on C(jti) with hi, < fi < /i c . Thus, the momentum 
distribution reaches a steady state as n — > oo and their 
mean energy (E) saturates to (E) . 

In Fig. Ufa), the broken curve (blue) shows (E) s for 
the classical system. As this figure shows, the mean en- 
ergy of the system increases with time and asymptotically 
approaches (E) . In the semiclassical regime, we expect 
a similar behaviour for the quantum average and this is 
shown as dashed curve in Fig. [SJa). The small differ- 
ence in saturated values of quantum and classical mean 
energies can be attributed to the finiteness of Planck's 
constant which makes its effect felt as e increases. 

Further, Fig. [8] (b,c,d) also shows the classical mo- 
mentum distribution f n (p) and its quantum analogue 
F n (p) = \ip{p,ri)\ 2 for the same set of parameters after 
evolving the system for n = 250, 275 and 300 kicking pe- 
riods. Probability distribution in position representation 
(not shown here) reveals that at n = 250, 275 and 300 the 
probability density in between the barriers is negligible. 
Nearly identical distributions in Fig. |5^b,c,d) mark the 
existence of steady state. Notice that small departures 
from semiclassical regime is also visible here in the form 
of slight difference between classical and quantum distri- 
butions. For the energy saturation effect, complete chaos 
between barriers is not essential. If some sticky islands 
are present between the barriers, the saturated classi- 
cal and quantum distributions as n — > oo will display 
a non-zero component in the finite well region. These 
non-chaotic components tend to remain localized and will 
never escape out. 
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FIG. 8: (Color Online) (a) Classical (solid line) and quantum 
(dashed line) mean energies as a function of time n for e = 
0.3. Other parameters are same as in Fig. [6] Numerically 
estimated value of {E) s for classical system is shown as a 
broken line. The triangles in the a;-axis are the times for which 
momentum distribution are shown in (b,c,d). Classical (solid 
line) and quantum (dashed line) momentum distributions at 
(b) n = 250, (c) n = 275 and (d) n = 300. Note that the 
distributions are nearly identical. 



C. Momentum Altering 

As demonstrated in section 5(B), when all the chaotic 
particles exit from the finite well region, a steady state is 
reached. One possible manifestation of this asymptotic 
state is the momentum filtering effect that occurs for cer- 
tain choices of parameters. It is possible to choose system 
parameters such that momentum distribution of escaped 
particles becomes narrow. Thus, any broad initial mo- 
mentum distribution at n = 0, after sufficient kicking 
periods, leads to a distinctly narrow momentum distri- 
bution. This is shown in Fig. |H1 In this figure, the initial 
conditions are uniformly distributed in the chaotic layer 
lying in between the barriers. This chaotic layer also en- 
sures that the final result is independent of the details of 
the initial distribution. The figure shows the momentum 
distributions fnoip) (classical) and i*Voo(p) (quantum) 
plotted for n = 700. By this time, a large fraction of 
particles have escaped from the well and the distribution 
has become bimodal with distinct peaks near —p c and 
p c . This shows that the double barrier structure, in pres- 
ence of the kicking field, acts as a momentum filter. We 
obtain filtering effect for a range of kick-strengths (not 
shown here) and observe that with decrease in e, the two 
bands in bimodal distribution become narrow. However, 




FIG. 9: (Color Online) (a) Classical (solid line) and quan- 
tum (dashed line) momentum distributions at n = 700 are 
displayed for V = 0.5, 6 = 0.2, R = 0.5, e = 0.1. For quan- 
tum simulation, h a = 0.0025. The initial distribution at time 
n — is a uniform distribution, the rectangular curve shown 
in blue. See text for details. 



the time at which system approaches steady state corre- 
sponding to this bimodal distribution becomes very large. 
Indeed, since it is experimentally possible to design bar- 
rier heights of desired choice, it will be be possible to 
use double barrier structure to produce filter with de- 
sired value of p c . From Fig. ® we note that the escaped 
particles follow extremely close set of invariant curves 
and their speed, averaged over time, will converge to the 
winding number of the orbits involved. Hence, the speed 
distribution will have peaks of infinitesimal width at \p c | . 
It is pertinent to note that a momentum filtering effect 
based on a very different mechanism has been studied by 
Monteiro et. al. in the context of a variant of kicked 
rotor model [22] . 

We emphasise that all the dynamical features discussed 
in section 5(A,B,C) can be explained on the basis of (i) 
co-existence of diffusive and non-diffusive regions which 
exists because the non-KAM nature of the system affects 
the dynamics differently in different phase space regions, 
and (ii) presence of KAM curves through out the phase 
space outside the double barrier region. Hence, all the 
dynamical features can be attributed to the interplay be- 
tween the KAM and non-KAM behaviour of the system. 



VI. DISCUSSIONS AND SUMMARY 

In summary, we have presented primarily numerical 
results of the dynamics of non-interacting particles in a 
double barrier structure acted upon by periodic kicking 
field. This model differs from the paradigmatic kicked 
rotor. This is essentially a non-KAM system and hence 
chaotic dynamics sets in for even for infinitesimal excur- 
sions from the integrable limit of kicking strength e = 0. 
Further, this displays non-equilibrium steady state and 
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classically induced suppression of energy growth in the 
semiclassical regime. This is in contrast with the classi- 
cal kicked rotor that displays diffusion only for e >> 1 
and its quantum version arrests this through dynamical 
localisation, an outcome of quantum interferences. 

Some of the earlier works on the double barrier type 
potential have considered it as a scattering problem, in a 
different setting with a drive term. For example, see ref- 
erences [23| . An incoming wavepacket hits the left barrier 
(see Fig. [lj and tunnels in to it and, depending upon the 
parameters chosen, some or all of it emerges out of the 
right barrier. This mechanism requires purely quantum 
effects such as tunneling and in this work we have deliber- 
ately avoided them to focus on the semiclassical regime. 
Since tunneling probability is nearly zero in this semi- 
classical setting, any initial distribution placed anywhere 
outside the barriers (\x\ > a + b) will continue to evolve 
on the KAM like invariant tori. However, based on the 
results obtained in this paper, we can speculate about 
the case when quantum effects come into play. Tunnel- 
ing will allow a wavepacket to enter the through the left 
barrier and non-KAM chaos will ensure that it gets dis- 
persed. But now, the wavepacket can tunnel out through 
the right barrier. This scenario could potentially lead to 
an interesting competition between above barrier cross- 
ings and tunneling. Another interesting case relates to 
periodic version of this model which can also be used for 
directed transport. We are pursuing these questions and 
will be reported elsewhere. 

The dynamical features in our model such as the non- 
equilibrium steady state and classically induced energy 
growth suppression are of current interest in the gen- 
eral context of transport and localisation especially for 
interacting systems such as the Bose-Einstein conden- 
sates. Recently there have been several experimental re- 
sults that point to classical features suppressing energy 
growth of condensates (Til ]. Typically, in such experi- 
ments, condensates are released from a confining poten- 
tial and their expansion in a disordered potential is stud- 
ied. When chemical potential fi < Vo, where Vo is the 
strength of disorder, condensates are classically reflected 
from the fluctuations of the disordered potential effec- 
tively localising the condensates. In our model, particles 
are neither interacting nor there is any disordered po- 
tential. However, the non-KAM chaotic dynamics and 
KAM like invariant curves provide the essential ingredi- 
ent for the suppression of diffusion. Even as the particles 
are transported in the position space their energy absorp- 
tion is restricted as t — > oo by KAM like structures. Such 
studies form an important background to understand and 
clearly distinguish similar quantum phenomena like the 
Anderson localisation from the classically induced ones 
and also to explore the connections between interactions, 
localisation and disorder. 

Quantum chaos in double barrier potentials have been 
studied before experimentally using GaAs/AlGaAs het- 
erostructures [24] though not with a periodic kicking 
field. In these experiments electrons tunnel through the 



double barrier potential and chaos is induced within the 
barriers due to the field created by the charge accumula- 
tion in the well 24]. Since resonant tunnelling plays an 
important role in this experiment, this can be regarded as 
being quantum in nature without classical analog. The 
double barrier system in Eq. [T] could be used with reso- 
nant tunnelling to study purely quantum effects as well 
though in the present work we have primarily explored 
the classical and semiclassical features. The foregoing ar- 
guments also imply that the system can also be realized 
experimentally in a laboratory. The cold atoms in op- 
tical lattices is the testing ground for variants of kicked 
rotor. An experimental set-up involving cold atoms, op- 
tical lattices with double barrier heterostructures should 
be possible. 

Currently there is considerable interest in the excit- 
ing field of chaotic ratchets [6j. Generally, ratchets are 
systems with broken spatio-temporal symmetries from 
which directed transport can be obtained even in the ab- 
sence of a net bias. There have been several proposals 
and at least one experimental realization for a chaotic 
ratchet in the last few years. The system presented in this 
work lacks the spatial periodicity required of a ratchet. 
But the kicking potential, being sinusoidal, is already 
spatially periodic. Further, from a theoretical perspec- 
tive, it is not difficult to have spatially periodic double 
barrier structures. Then, it might become possible to 
realize ratchet dynamics in this system. All the exist- 
ing chaotic ratchet proposals are based on systems that 
obey KAM theorem. The model presented in this work 
might lead to new ways to use non-KAM type dynamics 
for deterministic, directed transport. 



Acknowledgments 

The authors acknowledge many useful discussions with 
Dilip Angom during the course of this work. Numerical 
calculations for quantum system are carried out on PRL 
3TFLOP cluster computer. One of the authors (HP) 
thanks Manjunatha of CDAC and computer center staff 
for suggestions and assistance in efficiently using cluster 
computer. 



APPENDIX A 

Consider a particle that evolves on an invariant curve 
of the standard map C 5(115), approaches right barrier at 
x w = Rtt with p > p c during its motion after n th — kick, 
crosses it and exits on to another invariant curve of stan- 
dard map Ce(fj,e). In this appendix, we show that as the 
width of the barrier b — > 0, 65(^5) — > Cq(^). 

After the particle crosses the interface at x w and if At 
denotes the time it will take to cross the barrier region 
of width b, then At — > if b — > 0. Hence, the probability 
that a particle will experience the next kick while crossing 
the barrier will also tend to zero. Hence we can assume 
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that the particle does not experience a kick while crossing 
the barrier. In such a situation, the particle will face only 
two discontinuities between nth and (n+l)th kick. Thus, 
k = 2, B\ = x w and B2 — x w + b. From our assumptions, 

X J ) lie on C 5 (u 5 ), and ( X J ) will lie on C 6 (u 6 ). 

Pn J \ Pn J 



Pn 



Ul[ Pn 



fn 



Similarly, 



n 2 




Substituting for x\ and pi from Eq. IA1I in Eq. IA21 we 

get, 



? 1 - 



F 



Using b — > 0, we get 
C 5 (/x 5 ) -> C 6 (« 6 ) or u 6 - u 
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Pn 



(A3) 



This implies 



APPENDIX B 

We show that for certain special choices of (R, <fi) , re- 
flection from the walls of potential V sq takes a state from 



invariant curve C+ to its symmetric counterpart C_, 
where C+ and C_ are related through reflection sym- 
metry about (0,0). Let 



Rir- 
-Rtt 



In 

■ rrm 



!,m£ integer 



(Bl) 



Then, x r = In and — x\ = vnn. Let ^ J ^ e 011 
Reflection from the right boundary at x r will take it to 



Pn 



Pn 



2ln - x\- 



(B2) 



on the invariant curve C. The spatial periodicity of 2n 
in the standard map implies that 



2ln 



x'- 1 ) mod (2tt) 



~Pn 



(B3) 



is on C . Since 



-Pn 



is on C_ and CL is unique, 



we have C = C_ . Thus, the effect of reflection from the 
right boundary at x r is to take a state from C+ to C_ if 
Eq. IBll is satisfied. Similarly, the effect of reflection from 
left boundary at — xi is to take a state from C_ to C+. 
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